Sign-reversal of drag in bilayer systems with in-plane periodic potential modulation 
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We develop a theory for describing frictional drag in bilayer systems with in-plane periodic potential 
modulations, and use it to investigate the drag between bilayer systems in which one of the layers is 
modulated in one direction. At low temperatures, as the density of carriers in the modulated layer 
is changed, we show that the transresistivity component in the direction of modulation can change 
its sign. We also give a physical explanation for this behavior. 
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Throughout past decade there has been a great deal of 
experimental and theoretical activity in frictional drag 
in bilayer systems, following the seminal experiments by 
Gramila et a/.EJ These drag experiments involve a dou- 
ble quantum well system where the layers are individu- 
ally contacted by ingenious fabrication techniques. The 
barrier between the wells is made thick enough to sup- 
press tunneling but thin enough to allow significant in- 
terlayer interactions. An average current density ji is 
driven through layer 1 and circuit is kept open in layer 
2, so that J2 = 0. The interlayer interaction causes the 
electrons in layer 1 to drag along the electrons in layer 
2, and hence a counterbalancing electric field E2 forms 
in layer 2 to maintain a zero net j2- The transresistivity 
tensor *p*2i, defined by E2 = /?2iji, can be extracted ex- 
perimentally and can reveal important information about 
the properties of the effective interlayer interactions, the 
individual layers and the coupled bilayjer system. 

Since Gramila et aZ.'s original workd, which was done 
on a closely spaced electron-electron system at low tem- 
peratures without an applied magnetic field, many varia- 
tions on the theme of the original experiments have been 
performed. For instance, drag has been measured in 
electron-holds, and hole-holea systems, in widely sepa- 
rated layersffi—ajid in the presence of a perpendicular 
magnetic field. uu u Very recently, low-density systems 
have been studied to probe the suggested metal-insulator 
phase transition in stcoiigly correlated disordered two- 
dimensional systems.BEj Another modern trend— is to 
examine mesoscopic effects in Coulomb dragJlinlj. In 
general, drag without an applied i?-ficld is reasonably 
well understood within the frapewc-rk of a standard 
weak-interlayer coupling theory.liinllj The theory suc- 
cessfully accounts for several unusual features such as 
large enhancements in the transresistivity (up to an or- 
der the magnitude; some intriguing discrepancies ham- 
ever do persist for the most [dilute systems studiedEj) 
due to inttalayer correlationstH and plasmon mediated 
scattering. HLj On the other hand, the understanding of 
magnetodrag (i.e., drag in the presence of a perpendicu- 
lar £>-field) in bilayer systems is less complete, and sev- 
eral puzzling experimental results remain unexplained. 



For instance, under certain circumstances, the diagonal 
terms in the magnetotransresistivity (p%f and P21) nas 
been observed to reverse sign when the chemical poten- 
tial iSrGhanged in one layer while being kept fixed in the 
other .El This sign reversal with changing chemical poten- 
tial (which incidentally has not been observed at B = 0) 
cannot be obtained from magnetodrag calculations using 
the self-consistent Born appipximation,EjE3 and despite 
recent theoretical progress,E3 a fully satisfactory expla- 
nation of this phenomenon is not yet available. 

In this paper, we suggest that a reversal of the sign 
of the transresistivity is possible at B = in bilayer 
systems that have periodic potential modulations in the 
plane of the layers. The periodic potential modulation 
creates mini-bands, and the charge carriers can evolve 
from electron-like to hole-like behavior with a relatively 
small change in the density. Furthermore, for systems 
that are modulated in one direction, it is possible to ob- 
serve "skewed drag" (i.e., non-zero off-diagonal elements 
of P21), implying that the electric field response in the 
drag layer is in a different direction from that of the driv- 
ing current. This demonstrates the important role band- 
structure, plays in determining the transresistivity of the 
system.cJ We note that experiments on two-dimensional 
electron gases with strong potential modulations in one 
direction have already been reported in the literaturecJ, 
and hence we believe that the theory described below is 
amenable to experimental tests in near future. 

To investigate dragixL-these modulated systems, we use 
the Kubo formalismEj'EEl to calculate the transconductiv- 
ity tensor <72i, which is related to the transresistivity by 
P21 = — P22<?2i'pii in the weak interlayer coupling limit. 
In this method, the transconductivity is expressed as a 
current-current correlation function, which can be calcu- 
lated with standard perturbation theory techniques. 

The Hamiltonian of the system is H = Yli=i 2 + 
H12, where Hi is the Hamiltonian of layer i and 
H12 is the interlayer interaction term. We assume 
H21 is due to Coulomb interactions, so that H12 = 
A' 1 £^i(q)"2(-q)Vi2(q), where n 4 (q) and V 12 (q) are 
the Fourier transforms of the density operator and the 
interlayer Coulomb interaction, respectively. 
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We define, within the-ivlatsubara formalism, A to be 
the correlation functioned 

Jo Jo 

x (r T J(0)n(q,T)n(-q',r'))o. (1) 

For systems that have a periodic potential modulation 
with reciprocal vectors G, only q— q' = G terms are non- 
zero. Expanding in powers of V12, the first non- vanishing 
term for 021 in the dc limit is the second order term. We 
obtain 

<# = riE E Mq)M-q + GiM GlG2 

q G!G 2 

X / ^A^(q,q + G 2 ; W + i O+, W -zO+) 

x A7(q,q + G 1 ;w + zO + ,w-iO + )[-a w n B (w)] (2) 

Evaluation of A(q, q+G; ui+iO + , uj — i0 + ) is analogous 
to Ref. [h]. In this paper, we assume the system is in 
the weak scattering limit, allowing us to ignore vertex 
corrections at the charge vertices. Then, one obtains 

A(q, q + G, uj + i0+, uj - i0+) 

= -jYl [ V «'k+q r tr,n'( k + q) ~ V n kT t r,n (k)] 
krm' 

x [(n F (e nk ) - n F (e„ k - uj)] S (e„ k - e„'k+ q - u) 
xn(k + qn', kn; q) rj(kn, k + qn'; q — G). (3) 

Here, v„k is the band velocity, e is the energy, n F (e) = 
[exp(/3(e — (j,)) + (p, is the chemical potential), r tr 
is the transport time, r?(k' n', krt; q) = (k' n'\ exp(— iq • 
r)|kn) (n,n' are the band indices)]. The 0*21 obtained 
using Eq. (g) and the weak scattering result, Eq. (||), 
can also _he derived from the semiclassical Boltzmann 
equationO 

A complete calculation of drag, using (|^) in (||), is 
an arduous task, requiring a numerical evaluation of the 
band-structure(s) of the layers, calculation of the ma- 
trix elements ry(k'n', kn; q), and summation of different 
bands n, and reciprocal lattice vectors G. For incom- 
mensurate lattices one always has Gi = = G2, and in 
the remaining part of the paper we assume this to be the 
case. The other technical steps do not pose conceptual 
difficulties, and in the present context we find it appro- 
priate to consider simplified systems where to a certain 
extent analytic progress can be made, and for which the 
physics is transparent. 

The central issue of this paper is the sign reversal of 
the drag. We demonstrate this first for a 1-dimensional 
model, neglecting interband processes and the momen- 
tum dependence of the transport relaxation time. For 
this case, the correlation function (||) becomes A = 
T tr F(q,uj), where 



F(q,u) = — -^K- v k+q ) 
u k 

x [n F {e k ) - n F (e k+q )\ 8 (e k - e k+q - uj) 
= ^sign(w fei - v kz+q ) [n F (e ki ) - n F (e ki+q )} , (4) 

ki 

where ki are the solutions of e ki — £ k% + q — uj = 0. For 
illustrative purposes, we consider a cosine-band, e k = 
—h 2 /(ma 2 ) cos ka, for which there are two (or none) so- 
lutions, and one finds 

F(q, uj) = sign (v 2 - v{) [n F (-e 2 ) - n F (si) 

-n F (-Ei) - n F (e 2 )} ■ (5) 

At half-filling the chemical potential jj, vanishes, and 
making use of n FtlJ ,=o(— e) = 1 — 71jp,/*=o( e )j it is easy 
to see that the result (|^) vanishes identically. Thus, 
in an experiment where one of the subsystems is kept 
unchanged while in the other the chemical potential is 
moved through half-filling, the drag will change sign. 
While the above discussion is an important demonstra- 
tion of principle, it is necessary to also consider period- 
ically modulated two-dimensional electron gases, which, 
are the most commonly studied systems in this context .EZI 

A system which has an identical periodic modula- 
tion in both x- and y-directions is characterized by 
particle-hole symmetry, and it seems natural that the 
drag passes through zero when the two carrier species 
are matched. The experimentally most relevant sys- 
tems are those, however, where the modulation is only 
in one directiont3t3 (the strongest modulations have 
been achieved for these systems) and hence we choose 
the model system as follows: (1) There is a single 
band (the dispersion law and corresponding density of 
states are illustrated in Fig. 1) with a tight-binding dis- 
persion relation e(k x ,k y ) — h 2 [l — cos(k x a)]/ (m x a 2 ) + 
h 2 k 2 /(2m y ), and hence the velocity components are 
v x (k x ) = hsm(k x a)/(m x a) and v y {k y ) = fik y /m y . (2) r tr 
is k-independent. (3) The ipterlayer interaction ^12(9) is 
significant only for small qS3 

At low temperatures, it would appear permissible to 
expand in u), because dnsioj) / du in the integrand cuts 
off the higher uj contributions. Following this proce- 
dure yields an analytic expression for F. and the re- 
sulting drag resistivity obeys the familiar T 2 -dependeace 
known from unmodulated two-dimensional systemsEIE-X 
In this scheme, the transresisitivity diverges when // = 
2h 2 / (m x a 2 ) with the opposite sign from the low density 
pxx _ Thjs divergence, which is related to the divergence 
in the density of states at this energy (see Fig. 1), is 
unphysical because it only occurs in the experimentally 
unreachable T = limit. Nevertheless, it is interesting 
to observe that this approximation leads to a change of 
sign in the transconductance, and we should expect this 
behavior to be most prominent when [i = 2h 2 j(m x a 2 ). 

To cure this spurious divergence one must avoid expan- 
sions, and perform a numerical evaluation. As a starting 
point we have found it convenient to use 
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where 
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x [sin (k x - q x /2) a - sin (k x + q x /2) a] 

x [n F (k x - q x /2, k y0 (k x ) - q y /2) 
- n F (k x + q x /2,k y0 (k x ) +q y /2)} , (6) 
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transresistivity changes sign as the density is increased. 
The anisotropy of the transresistivity tensor implies that 
one should be able to see Hall drag. Experimentally, it 
may be possible to fabricate the system investigated here 
by by_overgrowing a pair of quantum weUs.over a cleaved 
edge^a or using lithographic techniques SB 
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While the Kubo formula gives the transconductivity, it 
is often most convenient to express the results in terms 
of transresisistivity (this is the object usually recorded in 
experiments) p^i, whose components are given by 
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and analogously for the yy-component. The transresisiv- 
ity tensor has the additional advantage that is does not 
involve the transport relaxation times for the individual 
layers, as long as these are momentum independent. The 
computed transresisitivity is shown in Fig. 2 for four dif- 
ferent temperatures. The most important feature is that 
the drag indeed changes sign; the effect is most promi- 
nent for low temperatures, and densities close to a fully- 
occupied band. In Figs. 3 and 4 we compare the xx- and 
yy-components of the computed transconductances (the 
ra-component was used in calculating the results of Fig. 
2). We observe that the sign change does not take place 
for Oil i nevertheless an interesting double-peak structure 
emerges. 

An analysis of the several assumptions made in our cal- 
culations is now in place. We have assumed that the sys- 
tem only has one band. Clearly, this assumption breaks 
down when the density so large that Fermi energy sig- 
nificantly exceeds 2ft 2 / (2m x a 2 ) , because the carriers will 
start to occupy higher bands. We also have assumed 
that temperature is low enough that the inelastic mean 
free path is much longer than period of the potential 
modulation, a. For finite £- lrL the system is roughly di- 
vided into coherent regions of order £f n . If i- m < a, the 
electrons do not coherently feel the periodic potential, 
and the drag characteristics will be given by an average 
of the drag over the density fluctuation caused by the po- 
tential modulation. Since the system as a whole acts like 
a (nearly) uniform system in this case, effects described 
in this paper will not be observable at temperatures for 
which £- m < a. 

To summarize, we have developed a theory for drag in 
bilayer systems where there is a periodic potential mod- 
ulation. We have calculated the drag for the case where 
there is potential modulation in one direction in one of 
the two layers. We find that at low temperatures the 
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FIG. 2. The calculated transresistivity pff, as a function 
of the density for four different temperatures. 
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FIG. 1. The density of states for a two-dimensional system 
with periodical modulation in one direction. The inset shows 
the constant energy surfaces for the dispersion relation used 
in this work. The energy is in units of fcgTo = 2h 2 /(m x a 2 ). 
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FIG. 3. The normalized transconductivity aff, as a func- 
tion of the density for the same temperatures as in Fig. 1. 
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FIG. 4. The normalized transconductivity a v ^j , as a func- 
tion of the density for four different temperatures. 
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